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ABSTRACT 

Visual and nfrared measurements from the five-channel AVHRR on the NOAA-7 
satellite are used operationally to derive sea surface temperatures. The multi- 
channel data are used to perform daytime and nighttime cloud-detection tests, and 
The several atmospheric-window channels in the thermal infrared »re used to 
correct for atmospheric attenuation. Monitoring of the sea surface temperatui e 
product with buoy data indicates s. ability in mean bias (<+0.1°C) and rms differ- 
ence (0.6-0.8°C.) with little variation by season or geographic area. Global 
napping enables the derivation of monthly mean isotherms, monthly and annua! 
changes, and anomaly patterns relative to climatology. Problems have been assoc- 
iated with noise in the 3.7/jm window channel, and with the injection of substantial 
volcanic aerosol into the stratosphere by the El Chi chon eruption. Multi-channel 
sea surface temperature charts are used in studies of such phenomena as equatorial 
long waves and the recent El Nino episode. 

1 . I NTRODUCT | ON 

Sea surface temperature fields on several spatial and temporal scales are of con- 
siderable interest to oceanographers, metec^o iog i sts, and climatologists. The 
irregular distribution and often variable quality of conventional intake temperatures 
temperatures from commercial ships, and the paucity of higher quality XBT and buoy 
observations, has stimulated interest in the repetitive and comprehensive coverage 
afforded by Earth-orbiting spacecraft. 

Rojtine processing of infrared ( IR) data from the NOAA series polar orbiting 
environmental satellites for global distributions of sea surface temperature (SST) 
began late in 1972 (Brower et al., 1976). Although the radiometers at that time 
performed their primary task of nighttime cloud imaging satisfactorily, the IR mea- 
surements were rather noisy and thus cloud-filtering was cumbersome and not as 
effective as required for extraction of SSTs. Furthermore, the single atmospheric 
"window" channel in the thermal-IR part of the spectrum did not permit adequate 
correction for atmospheric attenuation, chiefly by water vapor, which is * ry 
appreciable in tropical regions. Subsequent attempts to improve upon boti, cloud 
detec f ion and atmospheric corrections by incorporating data from atmospheric sounders 
aboard the spacecraft were only partially successful (Walton et al., 1976). Although 
still based on IR measurements from a single window channel, higher-resolution, lower 
noise- 1 eve I radiometric data from the Advanced Very High Resolution Radiometer (AVHRR) 
on the TIROS-N generation of NOAA operational satellites (Schwa lb, 1978), together 
with better sounder measurements, enabled further improvements in SSTs (Walton, 1980). 

Electrical interference in the 3.7jim window channel of the first AVHRR (on 
TIROS-N) hindered early development and implementation of multiple-window techniques, 
but simulations and preliminary testing with measurements from the second AVHRR (on 
N0AA-6) were quite encouraging (McClain, 1980). The third AVHRR (on NOAA-7) has 
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three IR window channels, and opera Mona 1 processing of multi-channel SST (MCSSD 
commenced in November of I 9b I (McClain et a I, , 1 9H 3 ) . 

This paper will give brief discussions of NOAA's MCSST processing, validation 
of MCbSTs, some global and regional products, and recent improvement activities. 
Use of MCSST analyses in atmospher i c/oceanograph i c research, and some gaps in our 
knowledge or understand i ng , will be addressed. 

2 . DESCRIPTION Oh MUL I I CHANNb L SSI (MUSS!) ML I HOD 

Spatial resolution of the AVHKR at nadir is l.l km locally and nominally 4 km 
regionally and globally. Measurements are made in the visual (0.58-0.68j;m) , re- 
flected IR (0.725-1 . lym) , and in three emitted-IR windows (3.55-3.93, 10.3-1 1.3, 
and I I .5- 1 2.5/im) . Each orbital swath is about 2500 km wide on the Earth and each 
geographic area is viewed near 0300 and 1500 local time. Only the IR channels 
have provision for onboard calibration. 

2.1 Cloud Detection Tests 


Bi-directional reflectance measurements from the visual or reflected-IR chan- 
nels are primary in the daytime cloud tests. Visible cloud threshold terts take 
advanfa e of the ocean/atmosphere reflectance being very low in the absence of 
clouds or glitter (specular ref I ectance) • These Thresholds have been determined 
as a function of solar zenith angle, satellite zenith angle, and azimuth of the 
viewed spot. When clouds partially fill a scanspot of the radiometer, reflectances 
from neighboring scanspots will generally differ by more than the low noise level 
of the measurements. This provides the basis for an additional test, esentially 
requiring spatial uniformity. Because both visible threshold and visible unifor- 
mity tests can erroneously fail in g I i tter-contami nated areas, the option exists 
to use the nighttime IR-un i formi ty test in the daytime as well. 

Nighttime cloud tests, of necessity, are based on data from the thermal-IR 
channels of the AVHRR. Threshold tes fs utilize the highest and lowest historical- 
ly observed ocean surface temperatures ; and a spalial uniformity test takes advan- 
tage of the high spatial resolution and low noise levels of the IR data taken in 
conjunction with the expectation of low horizontal gradients in both the SST and 
atmospheric attenuation. These assumptions are satisfied in most ocean areas and 
meteorological situations, but can fail at times, usually in coastal areas. Al- 
though these are the primary tests, a special test is used to detect low-altitude 
stratus clouds with uniform cloud top temperatures and that fill the radiometer's 
field of view. The reflectivity of opt i ca II y-th i ck water droplet clouds is sig- 
nificantly greater at 3.7 than at lljjm, and this has the effect of making the 
T 3 . 7*^11 difference much more negative than under any cloud-free conditions 
(McClain et al., 1983). Two other final filters can be used at night to remove a 
small residue of c , oud-confami nated cases. The first such test intercompares 
the MCSST values obtained from the dual-window, split-window, and triple-window 
equations (see next section), which have differing sensitivities to the presence 
of sub-resolution cloud contamination, and requires that they agree to within 
some limit, generally 0.5-1.0C. The second test requires that the triple-window 
MCSST not deviate markedly from the climatological SST (by more than 4-7°C, 
dependent on the climatological SST gradient). 

2.2 Correction for Atmospheric Attenuation 

Although the idea of using the different atmospheric transmittance character- 
istics of the several window regions in the infrared to devise an atmospheric at- 
tenuation correction scheme for satel I i te-der i ved SST goes back to at least Anding 
and Kauth (1970), the theoretical basis was explicated by Prabhakara et al. (1974) 
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and McMillin (1975), with the former being the first to verify the approach with 
radiometric measurements from space. A number of AVHRR simulations using model 
atmospheres and atmospheric transmittance models have been performed, among them 
Deschamps and Phulpin (1980), McClain (1981), and Barton (1983), and all found 
highly linear relations between the atmospheric correction and various combina- 
tions of channel brightness temperatures or brightness temperature differences, 
with standard errors of estimation generally falling in the 0.1-0.3°C range. 

NOAA's operational MCSST equations from NOAA-7 were derived from the following 
simulation equations ( °K in, °C out), which are based on a seasonally and geogra- 
ph ically diverse set of 64 atmospheric soundings from marine areas of the world 
and atmospheric transmittance models (Weinreb & Hill, 1980): 


4 

* 

4 , 


t 3/4 

- Til 

+ 

1.4887(T 3> 7-T n ) - 271.83 

(dua 1 -w i ndow) 

(D 

T 4/5 

= T 11 

+ 

2.4917(T n -T 12 ) - 273.48 

(sp! it-window) 

(2) 

T 3/4/5 = T 11 

+ 

0.95321 (T 3>7 -T ; 2 ) - 272.54 

(tr i p 1 e-w i ndow) 

(3) 


No explicit allowance for aerosol is incorporated in this set of equations, and the 
IR reflectance, the satellite zenith angle, and the sea-air temperature difference 
are a I I assumed to be zero. 
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3. VALIDATION OF MCSST RETRI EVELS 


t 


Initial testing of Eqn. (I) with N0AA-6 data (McClain, 1981) revealed temper- 
turedependent biases of up to 2 # C and a root mean square (rms) difference of 1.1°C. 
A relatively small set of expendable bathythermograph (XBT) and moored buoy 
measurements was employed to define temperature-dependent bias corrections for 
initial incorporation in the several MCSST operational algorithms (McClain et a I., 
1983). More extensive sets of drifting buoy observations have been collected 
and matched to satellite measurements over diverse geographic areas n d all seasons 
to determine the coefficients for the current operational equations v Strong and 
McC lain, 1 984 ) , g i ven below: 


t 4/5D = 1 *0346T , t 

+ 2.56(T n -T, 2 ) - 283.21 

(day) 

(4) 

T 3/4 = 1e0008T|j 

+ 1.50(T 3>7 -T 11 ) - 273.34 


(5) 

t 4/5N * 1 .0350T 1 1 

+ 2 • 58 < T 1 1 — T 1 2 ) - 283.18 

(night) 

( 6 ) 

T 3 / 4/5 = 1.0170T n 

+ 0.97(T 3 . 7 -Ti 2 > - 276.58 


(7) 


Validation of these equations with respect to independent satellite measurements 
and data from a totally different set of drifting buoys is given by the statistics 
in Table I. It should be noted that although MCSST equations involving measure- 
ments from the 3.7j»m channel can be used in the daytime, care must be taken to 
avoid sun glitter areas. Reflected solar radiation there elevates the brightness 
temperature above the value corresponding to emission alone. 
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Table 1. Validation of MCSSTs with respect to drifting buoy data 


Equation 


Sp I i t-w i ndow 
day n »ght 


Triple-window 
n iqht 


Dua I -w i ndow 


-0.02 -0.08 


- 0.01 


+ 0.02 


4. GLOBAL MAPPING OF MCSST 

Using satellite zenith angles of up to 45°, 12 hourly global coverage of 4-km 
AVHRR scanspot data is used to process hundreds of thousands of 11x11 scanspot 
arrays, called targets, over the oceans daily. Cloud tests and atmospheric 
corrections are applied to 2x2 unit arrays (8x8 km) within the targets, targets 
being centered about every 25 km. Using the most centrally located (day) or the 
warmest (night) of the cloudfree unit arrays within the target, 15,000-20,000 day- 
time and about 20,000 nighttime MCSST retrievals are made each day over the globe 
The less constrained daytime cloud tests allow observational densities of up to 
25 retrievals per target as an option. 

A variety of contoured MCSST fields is generated from the 8-km resolution 
retrievals, which are nominally spaced every 25 km or less depending on cloud dis- 
tribution and daytime observational density mode. All the daily (daytime and 
nighttime) MCSST retrievals are composited onto weekly 100-km (global) and 50-km 
(regional) grids, and a 14-km grid for coastal areas is under testing. The search 
distance for each grid point is an inverse function of the SST gradient, being a 
maximum of 200 km for the 100-km grid. Each MCSST retrieval within the search area 
is weighted, the maximum weight being applied to those within 1/3 of a grid inter- 
val, and the weight decreasing with the square of the distance from there. The 
distance-weighted average retrievals for each grid point each day are composited 
into a weighted average for the week by further weighting each retrieval as an 
inverse function of its age. 

The final weekly fields are contoured at 1°C intervals with the isotherms 
being dashed in areas where the grid-point MCSSTs have not been undated within the 
past seven days. Figure 1 is an example of a 50-km regional chart in the eastern 
equatorial area of the Pacific. 

All MCSST retrievals for the month are collected and averaged in 2.5° latitude/ 
longitude bins fcr purposes of generating fields of monthly means, departures from 
hi c+ orical means, (i.e., anomalies), month-to-month and year-to-year changes, etc. 
Figure 2 is a portion of the monthly anomaly chart for September 1982 during the 
early stages of a recent El Nino episode. 

5* MCSST IMPROVEMENTS 

Although the basic bias-corrected MCSST algorithms have been performing ade- 
quately under most conditions, certain situations are in need of further study. 

Among these are corrections for very large satellite zenith angles, for very 
large sea/air temperature differences, and for severe aerosol conditions. 
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Figure 1. 50-km MCSST analysis 
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Figure 2. Portion of global MC! 



5. I L.arqe Satellite Zenith Angles 
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Although the longer path length through water vapor at higher satellite zenith 
angles (sza) results in larger brightness temperature differences, and thus larger 
atmospheric corrections, it can be shown by simulation that this implicit compen- 
sation is inadequate when sza >^45°, particularly in the moistest atmospheres. 
Adding a term of the form K(Tj - TjMsec 0-1) where Tj and Tj are brightness temp- 
eratures in two of the window channels, 9 is the sza, and K is a constant, enables 
improved corrections at sza values of 45-60°. Table 2 below compares the satellite 
minus buoy/XBT statistics for 22 daytime (eastern North Atlantic) and eight night- 
time (Med i terranean ) matchups in March 1982 using Eqn. (2) and the following NOAA/ 
NESDIS splif-window simulation equation ( L I ewe I ! yn- Jones , 1983): 

T 4/5sza = T 11 + 2.346(T n - T 12 ) *■ 0.655(Tj 1 - T 12 Hsec 6 -1) -273.30 (8) 


ile 2. Sp 1 

it-window MCSSTs 
(Llewel 

with and without 
lyn-Jones, 1983) 

sza correction terms 

Equation 

N 

bias 

rmsd 

(2) 

28 

-0.44 

0.67 

(8) 

28 

-0.15 

0.44 


5.2 Large Sea/Air Temperature Differences 

It can be demonstrated by simulation that sea/air temperature differences 
(saAT) of to 5-!0°C, particularly if positive (sea warmer), result in very little 
error in retrieved MCSST when the algorithms based on saAT=0 are used. The very 
large negative saAT values, however, that are common in some areas (e.g., the Great 
Lakes and just off the NE coast of the USA in late spring and early summer) are 
associated with significant errors when certain of the MCSST equations are used. 

It is clear from the split-window and triple-window simulations in Fig. 3 that 
the latter equation is the preferred one in such situations, although the usual 
care to avoid specular reflection (glitter) areas must be taken in the daytime. 
Table 3 gives the bias and rms differences when the surface temperature is predict- 
ed by Eqn. (3) or Eqn. (4) using simulation data sets correspond i ng to T s -T a = -15, 
0, and + 1 5°C. 


Table 3. 

Simulations with 

MCSST equations 

assuming 

W0 



Sp 1 i t-w i ndow equat i on 

Triple 

-window equat i 

Ion 

Tc ~ T» 

N bias 

rmsd 

N 

bias 

rmsd 

-15 

67 +1.30 

1.86 

67 

+ 1.04 

1.10 

0 

67 +0.00 

0.38 

67 

-0.01 

0.07 

+ 15 

67 -0.63 

1.76 

67 

-0.05 

0.37 
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Figure 3. Split-window (T 4 / 5 ) 
and triple-window (T3/4/5) at- 
atmospheric corrections for sim- 
ulated extreme sea-air tempera- 
ture differences 
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5.3 Severe Aerosol Situations 


The dust veil from the eruption of £1 Chichon in April 1982 had a severe impact 
on MCSST processing for over a year, especially in the band from about 0-30°N lati- 
tude (Strong et al., 1983). Not only were daytime MCSST cloud tests failed because 
of reflected radiation from the volcano cloud, but nighttime retrievals were 
character i zed by substantial negative biases from increased attenuation by minute 
sulfuric acid droplets at very low temperatures. The Weinreb & Hill (1980) atmos- 
pheric transmittance models were used with a diverse set of 56 atmospheric sound- 
ings, with and without hypothesized El Chichon type aerosols, to determine simulat- 
ed brightness temperatures for the AVHRR window channels. The left two graphs of 
Fig. 4 show the simulated triple-window relation between the attenuation tempera- 
ture depression in Tjj and the temperature difference (Tj^-T^)* The aerosol-free 
data show the unique linear relation with extremely small scatter that is the basis 
of Eqn. (3) discussed in Sec. 2.2. The lower-left graph shows the same relation 
when two other aerosol concentrations are included, the 56 soundings being evenly 
divided into three aerosol classes. Each class exhibits a different relationship; 
i.e., no unique solution for SST appears. The graph at the right shows an alter- 
native triple-window relationship, the temperature depression in T 12 vs » Tj, 7 — T 1 1 ) , 
but with the same aerosol classes as at lower left. Although the scatter is great- 
er than in the aerosol-free case, there appears a unique relationship that is not 
sensitive to aerosol concentration. Preliminary testing of a volcano MCSST algori- 
thm based on this result was encouraging, but extensive testing was precluded by a 
severe noise problems in the 3.7yum channel (see Sec. 7). 

In addition, efforts are underway to extend the work of Griggs (1981) to severe 
aerosol conditions, and then to develop corrections to MCSSTs based on dynamic 
estimates of aerosol optical thickness from concurrent visual channel AVHRR data 
(Schwedfeger et al., 1983). Other extreme aerosol situations, such as Saharan 
dust outbreaks over the eastern subtropical Atlantic, also are known to affect 
certain MCSST retrievals, but investigations of these are just getting underway. 
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~ Figure 4 # Multi-window atmospheric corrections for simulated volcanic aerosols 


6. USE OF MCSSTs IN ATMOSPHERIC/OCEANIC RESEARCH 

Although available operationally for just two years, and plagued during 
much of that time with El Chichon effects and noise problems in the 3.7 m data, 
MCSSTs are seeing increased use both operat i ona I I y and in research. Legeckis* 
investigation of equatorial long waves in the eastern equatorial Pacific con- 
tinues to be greatly aided by the regional 50-km charts (Legeckis et al., 1983). 
Some of these westward-moving waves are readily apparent in the isotherms in Fig. 
1. The El Nino episode of 1982-83 was unusual in many respects (Rasmusson and 
Wallace, 1983), but among them was the intense positive SST anomaly that first 
began to develop along the equator between the Galapagos and the date line during 
the late summer of 1982. Few ships or other in situ observations we r e available 
in this area, especially during the early stages, and the MCSST analyses were the 
first operational charts to clearly depict the extent - and amplitude of the anomal- 
ously warm water in that region (Strong, 1983). The anomaly as it was in Septem- 
ber 1982 is seen in Fig. 2. MCSST techniques also were employed extensively by 
researchers involvd in the Gulf Stream Rings Project (Brown, 1983). 

1- CONCLUDING REMARKS 

"t' A' ‘'hough some of fhe physical processes are not really accounted for in the 

fairly simple theoretical models employed to date, especially those associated 
' with high concentrations of aerosols, the effects of most are sufficiently small 
and/or compensating that their net effect on the MCSST is also rather small, thus 
^ accounting for the stability of the good area-wide and season-wide accuracies 

measu r ed relative to drifting buoys. Similar matchup statistics have been obtain- 
ed daily on a global basis as part of the MCSST operational processing with re- 
spect to intake temperatures from so-called ships of opportunity. The statistics 
against ships tend to be poorer and more variable than those against drifting 
huoys (Strong and McClain, 1984), but this follows from the fact that ship SSTs 
matched with other ship SSTs (within 24 hours and 100 km) consistently have been 
found to have differences with a standard deviation of 1.5-2.0*C. There is the 
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possibility of increasing the accuracy of mu t t i -w i ndow bbTs yet further by adding 
further by adding multiple viewing-angle capability. The results of simulations 
at the Rutherford Appleton Laboratory in England have been encouraging (Llwellyn- 
Jones, 1983), and this approach is incorporated in the design of the Along Track 
Scanning Radiometer proposed for the ERS-1 satellite. 

A often asked question is the difference to be expected between the "skin 
temperature** sensed by a radiometer in space and a "bulk temperature** measured 
by an in situ sensor on a buoy (depth of 1 meter) or in a ships water- intake 
manifold (depth of perhaps 10 meters). This skin minus bulk temperature differ- 
ence should be 0.2°C or less with well-mixed upper layers; but in the daytime 
under a high sun and very light winds, the skin temperature can easily be elevat- 
ed several degrees above that at depth. 

A recurring engineering problem has been electrical interference in the 3.7^ 
m channel, the noise level increasing with time following launch of the space- 
craft. The less severe noise levels can be compensated for in the MCSST process- 
ing (e.g., using larger unit arrays) at the consequent expense of lower observa- 
tional densities and somewhat degraded accuracies from poorer cloud filtering. 
NOAA-7 noise levels had increased so much by June 1981 that it necessitated 
abandonment of this channel because of serious nighttime SST product degradation. 
Outgassing procedures were carried out on NOAA-7 for a week in September 1983, 
and this reduced the noise to very low levels again. Monitoring is being main- 
tained to ascertain whether the noise returns and at what rate of increase. 


Although serious attempts to map sea surface temperature globally from space . 

have been underway since early in the 1970*s, c n ly since late in 1981 when multi- t 

channel sea surface temperatures began to be processed operationally from AVHRR J 

measurements, have the coverage and accuracy achieved levels where oceanographers 
and meteorolog i sts alike are beginning to use sate M i te-der i ved SSTs in their ] 

activities. With good coverage provided by sate I I i te-based sea surface tempera- \ 

tures in the Southern Hemisphere oceans and other data-sparse areas, a truly 
global time series becomes feasible. j 
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